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The dynamics of transverse Neel domain wall in a ferromagnetic nanostrip in the 
presence of driving field, current and transverse magnetic field is investigated by the 
Landau-Lifshitz-Gilbert(LLG) equation with the adiabatic and non-adiabatic spin- 
transfer torques both analytically and numerically. The analytical expressions for the 
velocity, width, excitation angle and displacement for the domain wall are obtained 
by using small angle approximation along with Walkers trial function. The results 
show that the initial velocity of the domain wall can be controlled by the adiabatic 
spin-transfer torque and the saturated velocity can be controlled by the non-adiabatic 
spin-transfer torque and driving held. The large increase in the saturated velocity of 
the domain wall driven by current and held due to the transverse magnetic held is 
identihed through the presence of driving held. There is no impact in the saturated 
velocity of the domain wall driven by current from the transverse magnetic held. For 
the domain wall driven by the current in the presence of the transverse magnetic 
held, the saturated velocity remains constant. The transverse magnetic held along 
with current and driving held is more advantageous that the transverse magnetic 
held along with current for increasing the saturated velocity of the domain wall. The 
numerical results showed that the saturated velocity is increased by the transverse 
magnetic held with the irrespective of the directions of the driving held and current 
further it is higher and lower when the directions of driving held and current are 
antiparallel and parallel respectively. The obtained analytical solutions are closely 
coincided with the computed numerical results. 
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I. INTRODUCTION 

Domain wall motion in magnetic nanostructures driven by a magnetic field and cur¬ 
rent have variety of potential and technological applications including logic devices^, atom 
trapping^ and memory storage^”®. Both the field induced and current induced domain 
wall dynamics have been extensively studied experimentally^’^“^^ and numericallyThe 
transverse domain wall has higher velocity and less complexity in motion in the presence of 
magnetic field^^. In the presence of current or field the transverse Neel wall in nanostrip can 
be moved along the length axis and the average speed of the wall is almost linear with field 
or current below the critical value, the so called Walker limit of field or current respectively. 
The domain wall moves rigidly with slight excitation and distortion below the Walker limit 
whereas above the Walker limit, the motion of the domain wall is no more linear and the 
average speed of the wall reduces suddenly. These have been stuided analytically^*^’^^’^^, 
numerically^®’^^, experimentally^^and also through micromagnetic simulation^^’^®’^'^’^^“^°. 
The drastic decrease in the speed of the domain wall is due to the oscillatory behaviour that 
occurs in the structure of the wall, which depends on the cross sectional area of the strip. 
When the width and thickness of the strip are large, the domain wall oscillates between 
transverse'^^’’^"^ and vortex type domain wall, whereas in the case of the width and thickness 
are small, the transverse wall rotates about the length axis^^’^^’^^’^®. There are few meth¬ 
ods available to increase the speed of the domain wall and the Walker limit: introducing 
roughness in the stripmaking nanostrip with honey comb structure^^, applying oscilla¬ 
tory magnetic field^ ' , inclusion of perpendicular anisotropy under layer’^® and by applying a 
transverse magnetic field. 

The recent experimentaR®”'^^ and micromagnetic simulation^^“'^^ studies on field induced 
domain wall dynamics show that when a transverse magnetic field is applied parallel to the 
magnetic moments of the domain wall and it maintains the regular motion of the domain 
wall even for large applied field. This implies that, the Walker limit is increased and it 
corresponds to increase the speed of the domain wall. Also, the presence of the transverse 
magnetic field, increases the width of the domain walh^^a making asymmetry and 
twisting in the domain wall"^”^. However, when a transverse magnetic field is applied in a 
direction antiparallel to the direction of magnetic moments of the domain wall, the speed 
and width of the domain wall as well as the Walker limit decreases"^^’'^^. 
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Analytical and micromagnetic simulation study of current induced Bloch wall^^ motion 
in perpendicularly magnetized nanostrip and micromagnetic simulation study of Neel wall^® 
motion in parallely magnetized nanostrip show that there is an increase in the Walker limit 
as well as the speed of the domain wall. Inspite of the above developments in the case of 
field induced domain wall motion in the presence of a transverse magnetic field, the current 
driven domain wall dynamics along with transverse magnetic field is yet to be understood 
properly, especially for the transverse Neel domain wall. A systematic analytical study on 
the dynamics of Neel domain wall driven by current and field in the presence of a transverse 
magnetic held is not available in the literature. 

Motivated by the above, in the present paper, the dynamics of Neel domain wall in a 
ferromagnetic nanostrip driven by both the current and held in the presence of a transverse 
magnetic held are extensively studied through analytically and numerically. The analytical 
expression for the excitation angle of the magnetization inside the domain wall, the velocity, 
the width and the displacement of the domain wall are obtained. The paper is organised 
as follows: In Section II, we present the model and the governing equation of motion for 
the dynamics of domain wall in a ferromagnetic nanostrip. In Section III, the dynamical 
equation is analytically solved using Walker’s trial function and the dynamical parameters 
such as excitation angle, velocity, width and displacement are derived. The numerical results 
are obtained by solving the dynamical equation using Runge-Kutta-4 method and the results 
are compared with the analytical results in Section IV and the effects of current, driving field 
and transverse magnetic field on the dynamics of domain wall are discussed. The influence 
of the transverse magnetic field in saturated velocity is also observed. Finally, the results 
are concluded in Section V. 


II. DYNAMICAL MODEL FOR NEEL TYPE DOMAIN WALL MOTION 

Consider an infinitely long anisotropic ferromagnetic nanostrip with an easy axis of mag¬ 
netization along x-direction as shown in FIG.l. In the nanostrip, the uniformly magnetized 
left domain along -|-x direction and the right domains along -x direction are separated by a 
domain wall, where the magnetization can be represented by the vector M(a;, t). In FIG. 1(a), 
and J represent the driving field and current density applied along x-direction. The trans¬ 
verse magnetic field is applied along the positive y-direction. The arrows indicate the 
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(a) (b) 


FIG. 1. (a) A schematic representation of ferromagnetic nanostrip which is taken as our model 

with its easy axis of magnetization along x-direction. The current density J and external magnetic 
field Hrf can be applied in both the directions of x-axis and the velocity of domain wall depends 
upon the directions of field and current. The transverse magnetic field is applied along positive 
y-direction. 


direction of magnetization along the nanostip. In FIG.1(b), e^,, By and represent the unit 
vectors along x, y and z-directions respectively. It is assumed that the magnetization of 
the nanostrip varies only along the x-direction and the initial prohle of the magnetization 
of the model appears as shown in FIG.1(a). The angles 6{x,t) and ^{x,t) represent the 
orientation of magnetization with reference to the positive x-direction and its projection 
in the yz-plane making with the positive y-direction respectively. Physically the angle 6 
represents the spatial variation of the direction of the magnetization along x-direction and 
<h represents the out-of-plane excitation of magnetization vector in the strip. The Landau- 
Lifshitz-Gilbert(LLG) equation that governs the dynamics of the magnetization present in 
the strip in the presence of an externally applied fields and current is written as'^^ 


(9M(x, t) 

m 


a 


aM 


= -yMxHe// + —Mx^-—Mx (Mx 


dt Mr 


aM 

dx 


c 


( 1 ) 


M = (A4, My, M,); |M|2 = M^ + M^ + Mj = Ml 


Here, M(a:,t) represents the magnetization vector, 7 is the gyromagnetic ratio, a is the 
Gilbert damping parameter, Ms{= |M|) is the saturation magnetization of the magnetic 
strip and b = PJfiB/eMg, c = represent the magnitude of adiabatic and non-adiabatic 
spin-transfer torques respectively. Where, P is the polarization, J is the magnitude of 
the current density, ps is the Bohr magneton, e is the charge of the electron and ^ is 
the non-adiabaticity factor. In Eq.(l), He// represents the effective field due to different 
magnetic contributions including exchange interaction, magnetocrystalline anisotropy, the 


5 














R. Arun, P. Sabareesan and M. Daniel 


driving field, the transverse magnetic field and the demagnetization field. The first term 
in the right hand side of Eq.(l) represents the precession of the magnetization about the 
effective field He//, that determines the precessional frequency and conserves the magnetic 
energy. The second term, supports the damping of the magnetization due to dissipation of 
energy that takes place within the material. The third term, which represents the adiabatic 
spin-transfer torque corresponding to the reaction torque on the magnetization produced by 
the spatial variation of the spin current density^®. The last term represents the non-adiabatic 
spin-transfer torque, which corresponds to the reaction torque on magnetization due to the 
continuous space variation of spatially mistraking spins between conduction electrons and 
local magnetization. Adding all the above fields, the total effective field can be written as 



where A represents the exchange interaction coefficient and Hk represents the magnetocrys- 
talline anisotropy coefficient. H^, Hy and are the magnitudes of driving field, trans¬ 

verse magnetic field and demagnetization field respectively. 

The dynamics of magnetization in the domain wall can be understood by solving Eq.(l) 
after substituting Eq.(2) in it. While solving the dynamical Eq.(l), we also evaluate the 
dynamical quantities of the domain wall namely velocity, width, excitation of the wall from 
the plane of the strip and the displacement analytically. 

III. ANALYTICAL SOLUTIONS FOR THE DOMAIN WALL 
PARAMETERS 

As Eq.(l) is a highly nontrivial vector nonlinear evolution equation, it may be difficult 
to solve the same in its present form. Hence, we rewrite Eq.(l) in terms of the angles 9{x, t) 
and <h(x,t) by defining the component of magnetization in the polar form as follows. 


Mx = Ms cos9{x, t), 

My = Ms sin 9{x, t) cos <h(x, t), 

M 2 = MsSm9{x,t)sm^{x,t)- 


(3a) 

(3b) 

(3c) 
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On substituting Eqs.(3) in Eq.(l) along with c = ^b, we obtain the following set of equations 
for 9{x,t) and <h(x, t). 


d9{x,t) . d^{x,t) 

^ ^ +asm9{x,t)—= 7 


dt 


dt 


2A 

w. 


.d9(x,t)d^ix,t) m sd‘^^(x,t) 
2 cos 9{x,t) +sm9{x,t)- 


dx dx 


— Hy sin $(a:, t) — 2 ttMs sin 9{x, t) sin 2$(a:, t)] + b 


d9{x, t) 
dx 


+ ,^fesin 6 *(a;, t) 


dx"^ 
5<h(x, t) 
dx 
(4a) 


d9{x,t) . .d^{x,t) 

a—- sm9{x,t) - 57 ^= 7 


dt 


dt 


^ I d^9(x,t) 

Mo \ dx^ 


sin 9(x, t) cos 9(x, t) 


5$(x, t) 

dx 


— Hd sin 9{x, t) + Hy cos $(a:, t) cos 9[x, ~ 2 sin^ <h(a;, t)] sin 29{x, t) 


^d9{x,t) . n/ sd^{x,t) 

^ -sin 0 (a:,t)' ^ ^ 


dx 


dx 


(4b) 


Eqs.(4) describe the dynamics of domain wall in a ferromagnetic nanostrip in the presence 
of current, driving field and transverse magnetic field in terms of 9 and <h. The problem 
now boils down to solving Eqs.(4) for 9, $ in order to find the excitation angle, velocity and 
width of the domain wall. The same form of Eqs.(4) with and without transverse magnetic 
held in the absence of current has been studied and solved by Schryer et and Lu et 
respectively using the trial functions obtained from the static prohles 9{x/W) and $(a;), 
which are the spatial variations of 9 and <h when the domain wall is at rest and W is the 
width of the domain wall. Similarly, Li and Zhang^^ have studied the dynamics of transverse 
Neel wall in the presence of driving held and current, without transverse magnetic held, using 
the trial functions obtained by Schryer and Walker^®. In the present paper, the dynamics 
of transverse Neel wall driven by driving held and current in the presence of transverse 
magnetic held is studied by solving the Eqs.(4) using Schryer and Walker’s trial functions, 
which are given by 


9{x, t) = 9 


X 




W(t) J 

<h(a;, t) = <h(a;) + (j){t) U 


( x-X{t) \ 

V W{t) )■ 


(5a) 

(5b) 
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Where X (t) is the position of the center of the domain wall and U is a. step function which 
can be defined as 


U = 1, when, 



(6a) 


and 


17 = 0, when, 



(6b) 


The new assigned function (j){t) can be defined as the angle of the magnetization at the 
center of the domain wall, can be called as excitation angle. In the forthcoming sections, 
the trial functions will be obtained after deriving the static profiles of 6 and 

A. Static profile for the domain wall in the presence of transverse magnetic 


field 


Following Schryer and Walker’s^® procedure, it is required to hnd out the trial functions 
of 6 and $ from their static profiles and these trial functions are to be substituted in Eqs.(4a) 


and (4b) to determine the excitation angle, the velocity and the width of the domain wall. 
In the present work, the static profiles for 6 and $ are derived after the transverse magnetic 
held is applied when the current and driving held are switched oh. When a transverse 
magnetic held Hy is applied along the positive y-direction, it exerts a torque on the magnetic 


moments in the strip and changes their direction from the initial equilibrium direction to 


the new equilibrium direction in the duration of few piccoseconds and as a result of this, 
the magnetization inside the domains are tilted towards the positive y-direction from the 


direction of easy axis. Therefore, the value of 6 in the entire strip except at the centre of the 


domain wall changes to the new orientation in the presence of a transverse magnetic held 
and correspondingly the width of the domain wall increases. Whereas there is no variation 
in $ and it takes the value of zero for the entire strip even after the transverse magnetic 
held is switched on, because the transverse magnetic held is applied along the same plane 
in which the magnetic moments are present. Eventhough, the direction of magnetization in 
the domains changes, there is no spatial variation in the direction of magnetization inside 
the domains and thus the domain wall is static in the presence of a transverse magnetic 
held. Hence, the static prohles can be found in the presence of a transverse magnetic held 
and absence of current and driving held. In the absence of any external held and current. 
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the value of 6^ is 0 and vr in the left and right domains respectively and it varies from 0 to 
TT inside the domain wall from left to right. When the transverse magnetic field is switched 
on, the values of 6 in the left and right domains are changed and they can be represented 
as 9D and ti — 9d respectively. 


H 

A' 








► e 


FIG. 2. The spatial variation of magnetization vector of the nanostrip in the presence of a transverse 
magnetic field Hy applied along the positive y-direction and in the absence of a driving field and 
current. The angle 9 varies from 9£, to n —9£> from left edge to the right edge along the x-direction 
and the angle <I> is zero in the entire strip. 

Static profile 9{x) can be computed in the presence of transverse magnetic field once 
the current and driving field are switched off and the following equation is obtained after 
substituting = b = 0 and <h = 0 in Eq.(4a). 


£9 HkMs . ^ ^ , HyM, ^ ^ 

—t; -:— sm 9 cos 9 H-— cos 6* = U. 

dx^ 2A 2A 

After multiplying Eq.(7) by it is integrated with respect to x as 

fd9Y H,M, , 2 , , HyM, , ^ ^ 

I — I-— sm 9 H-— sm 9 = Ci, 


(7) 


\dx J 


2A 


A 


( 8 ) 


where Ci is the constant of integration and it is obtained after substituting ^ = 0 and 


9 = 9£)(ior the left domain) in Eq.(8) as 


= ^^sm9n -^^sm^9 


D- 


(9) 


A 2A 

The value is found from Eq.(4b) after substituting b,Hd = 0 and 9 = 9d,^ = 0 
corresponding to the left domain as follows: 


sin So = = sin-‘ . (10) 

9d from Eq.(lO) gives the angle between the direction of magnetization in the left domain 
and positive x-direction in the presence of transverse magnetic field. Knowing 9j), the angle 


9 









R. Arun, P. Sabareesan and M. Daniel 


between the magnetization and positive x-direction in the right domain is calculated as 
(vr — 6*£)). From Eq.(lO) one can understand that when Hy = Hk, 9 d becomes 7r/2, which 
means that all the magnetic moments will be aligned parallel to the direction of the applied 
transverse magnetic held, which leads to the disappearence of the domain wall. Therefore 
one can set the condition for the transverse magnetic held as Hy < Hk- By substituting 
Eq.(lO) in Eq.(9), Ci is obtained as 


Cl 


MsHl 

2AHk 


( 11 ) 


On substituting the value of Ci from Eq.(ll) in Eq.(8) one obtains 


dO 


= ± 


HkM, 


sin6' 


Hu 


( 12 ) 


dx V 2^4 

The positive sign in Eq.(12) represents the variation of 9 from 0 to vr in the strip along the 
positive x-direction(head-to-head domain wall) whereas the negative sign is corresponding 
to the variation from tt to O(tail-to-tail domain wall). In the present model (FIG.l), the 
head-to-head domain wall conhguration is considered to study the domain wall dynamics, 
and therefore in Eq.(12), only the positive sign survives. On integrating Eq.(12) we obtain 


: In 


Jtan|-Wl 


H, 


Hk ) 


1 - 


1- ftlanf + 


1 - 


HkM, 

2A 


X C*2, 


(13) 


where C 2 is the constant of integration which can be determined by substituting 9 = 7r/2 
and a: = 0 in Eq.(13). The result reads 


02 = 


In 


1 _ Ejl 

Hk 


1 + 

^ Hk 


1 - 


i-fe+vi+a 


On substituting O 2 in Eq.(13), one can derive the static prohle for 9 as 

ai -F 02 exp (I 7 ) 


9{x) = 2 tan 


-1 


(14) 


(15) 


02-^01 exp(|7) 

where, oi = + 02 = y^l + §+-\/l " ^ = ^o/y^l- (Sf) and Wc 


2A 


HkM'- ^ and Wo are the width of the domain wall in the presence and absence of trans¬ 
verse magnetic held respectively when current and driving held are switched oh. The Eq.(15) 
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represents the spatial variation of the angle between M and positive x-direction after the 
transverse magnetic field is applied while the driving field and current are switched off. In 
the absence of transverse magnetic field {Hy = 0), Eq.(15) reduces to 2tan“^(exp(x/hFo ))5 
which agrees with the result obtained by Li and Zhang'^^. The components of the mag¬ 
netization in the strip in the absence of current and driving field and in the presence of 
transverse magnetic held are obtained by substituting the Eq.(15) and $ = 0 in Eqs.(3). 


My 


Ms cos 


^2 tan ^ 

ai + 02 exp 


02 -b oi exp [^) 

) 

^2 tan“^ 

oi -b 02 exp {^) 

] 

02 -b Oi exp (^7) 

) 


(16a) 

(16b) 

(16c) 



X (nm) 



FIG. 3. (color online), (a) The spatial variation of 9 for Hy = 0, 200,400, 500 Oe in the absence of 
driving field and current, (b) The spatial variation of the components of magnetization {mx,my 
and niz) for Hy = 0, 200 Oe. 


The static prohle 6{x) and the normalised components of the magnetization m(= M/M^) 
is given in Eqs.(15) and (16) respectively have been plotted in FIG. 3 for various values of 
transverse magnetic held. When the strength of the transverse magnetic held is increased 
from 0 to 500 Oe, the value of 6 in the left domain and right domain increases and decreases 
towards 7r/2 respectively as shown in the FIG.3(a). The slope of 6 at the center of the 
domain wall decreases while increasing the Hy from 0 to 500 Oe. This implicitly refers to the 
expansion of the domain wall width. The spatial variation of the normalised components 


11 







































R. Arun, P. Sabareesan and M. Daniel 


of the magnetization m = M/M^ is shown in FIG.3(b). While the value of is zero 
everywhere, and My are varied spatially. The broadening of the hump in the plot 
corresponding to My due to the increase in Hy from 0 Oe to 200 Oe and it confirms that, 
the increase in width of the domain wall. 

B. Trial function for the moving domain wall 

The static profiles for <F and 9 after applying transverse magnetic field, in the absence 
of current and driving held, are given by <F(a:) = 0 and 0(x)(Eq.(15)) respectively. Now we 
compute the trial functions of $ and 9 corresponding to the moving domain wall driven by 
the driving held and the current using the earlier static conhguration. The applied current 
and driving held excert the spin-transfer torque and torque to the magnetic moments in 
the domain wall respectively and they create an energy diherence between the two domains. 
In order to reduce the total energy of the nanostrip the domain wall moves. Though the 
presence of driving held can change the orientation angle 9 in the left and right domains 
from the previous orientation angle and — 9d) respectively and this variation can be 
neglected when the strength of the driving held is small. The another orientation angle 
<h inside the two domains remains unchanged because the magnetic moments excite from 
the xy-plane and settle back to the same plane immediately due to damping. In the case of 
current along with transverse magnetic held, the orientation angels 9 and <h remains constant 
due to the absence of both the adiabatic and non-adiabatic spin-transfer torques within the 
domains. However the orientation angels 9 and <F varies inside the domain wall due to the 
torque and damping arise by the current and driving held. By neglecting the variation of 9 
and <F inside the domains due to driving held and current, the trial functions for 9 and <F 
corresponding to the moving domain wall can be obtained by substituting the Eq.(15) and 
<h(a;) = 0 in Eq.(5a) and Eq.(5b) respectively and the results read 



(17b) 


(17a) 
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The trial functions Eqs.(17) represent the modelled solution for 6 and $ corresponding to 
the moving domain wall driven by current and driving field in the presence of transverse 
magnetic held. Using the above trial functions the unknown quantities such as excitation 
angle((/)), velocity(n), width(lU) and displacement(X) of the domain wall will be derived in 
the forthcoming section. 


C. The dynamical quantities of the domain wall 


On substituting || from Eq.(4b) in Eq.(4a), we obtain 


d4) 

'm 


(1 + sin 6* I M, 


2A7, ^2$ aed^ d‘^e 

a sm 0 TTvr + 2a cos 0 - 7 —— - 7-^ + sm 0 cos 0 


dx"^ 


dx dx dx"^ 


7r 


—2a77rMs sin 2<h sin ^ ^ [Hk + 4:TtMs sin^ $] sin 26 + 'jHd sin 9 

(9$ 06 

-a'yHy sin $ — ^Hy cos <h cos 6 + {1 -\- a^)h sin 6 ——h (a — — 

(yJb (JJb 



2 


( 18 ) 


Similarly, by substituting ^ from Eq.(4a) in (4b), we get 


06 


2 A 7 , ^ 2 $ 06 0 <!> 0 ‘^e 

—2'ynMs sin 2<h sin 6 — ^ [Hk + 47rMs sin^ $] sin 26 — a'yHd sin 6 


Ox J 


5 $ do 

-jHy sin ^ + a'yHy cos ^ cos 6 — {a — ^)b sin 6* ——h (1 + a^)6 — [■ . 

(y Jb (J Jb 


( 19 ) 
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From Eqs.(17a) and (17b), one can derive the following identities at x = X{t). 


9{X,t) = n, 

de{x,t) 1 


'Hk-Hy 


W{t) V Hk + Hy^ 


= 0 , 


Hk 


IHk-Hy 


W{ty {Hk + Hy) ]j Hk + Hy 
1 fdX\ iHk - Hy 


dx 

d^e{X, t) 
dx"^ 

d^e{X,t) 

dx^ 

de{x,t) _ 

m “ ~ W{t) \dt J y Hk + Hy' 

d^{X,t) d(j){t) 

dt dt 

d^{X,t) _ ^ 

dx dx'^ dxdt 

After reducing Eq.(18) at x = X{t) and substituting Eqs.(20a)-(20h) in it, we get 


2 .d(j){t) 


(1 + a ) = 'y{Hd - a[27rMs sin20 + Hysincj)]) + 


{a-Ob Hk-Hy 


dt -r • \j Hk + Hy 

Similarly, the velocity of the domain wall is obtained from Eq.(19) as 

1 + 


v{t) = 


dX ■yW (t) lHk + Hy 


dt (1 + aO V Hk-Hy 


[27rMs sin 20 + aHd + Hy sin 0] — 6 


1 + 


(20a) 

(20b) 

(20c) 

(20d) 

(20e) 

(20f) 

(20g) 

(20h) 


( 21 ) 


( 22 ) 


The width of the domain wall is obtained after differentiating Eq.(18) with respect to x and 
reducing the differentiated equation at x = X{t) using Eqs.(20a)-(20h). 


W{t) = IFoi 


Hk 


Hy + Hk 


1 + 


AttM, 


S • 2 i 

- sm 0 —^ cos ( 

Hk Hk 


H. 


(23) 


and the ratio of the width to the initial width namely the width ratio of the domain wall is 
given by. 


Wjt) _ 

fR(0) V Hk 


AnMs . ^ Hy ■ 

1 H-—— sm 0 — -jO cos 0 

Hk Hk 


(24) 


Eqs. (21), (22) and (23) represent the excitation angle(0), velocity(x) and width(hF) of the 
domain wall driven by driving held and current in the presence of a transverse magnetic 
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field. From the equations (22) and (23), one can understand that the velocity and the width 
of the domain wall are the function of excitation angle. Hence, we need to solve the Eq.(21) 
in order to understand the domain wall motion. The Eq.(21) is a highly nontrivial nonlinear 
evolution equation and solving it analytically in the present form is extremely difficult. But it 
is possible to find out the exact analytical solution of Eq.(21) only when ^ 0,b = Hy = 0 
and Hy 0,h = Hd = 0 and for the other cases, the Eq.(21) can be solved using small angle 
approximation technique which would be discussed in the forthcoming sections. 

D. Exact analytical solution for the domain wall motion 





FIG. 4. (color online). The variation of (a) velocity, (b) width ratio and (c) excitation angle of the 
domain wall with respect to time driven by different driving fields below Walker limit(181.6 Oe). 
The corresponding inset figures show the variation above Walker limit. 


1. In the presence of driving field alone(Hd 0 and b = Hy = 0): 

To study the dynamics of the domain wall in the presence of driving field and absence of 
current and transverse magnetic field, Eq.(21) is rewritten as 

(1 + a^) ^ = 'y{Hd — 27rQM^ sin 20). (25) 

The solution 0(t) can be obtained by integrating the Eq.(25), is given as 



d(f) 

Hd — 27raMs sin 20 



7 

1 + 


dt. 


( 26 ) 
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For the initial condition is 0(0) = 0, the excitation angle 0(t) is solved from Eq.(26) as 
0(t) = tan“^ 


{G -|- 27iaMs) exp ■ 

!- 

In 

^TTOil^s ~\~G 

27raMs—G 

+ !&*)] 

>+G 

— 2'KaMs 1 

Hd (^exp < 

[- 


2 TT Q; iVfs -|-(j 

2naMs—G 


00 

i 


(pit) = tan’ 


when 

j 27raM, /, f 27raMs 


\ 


+ 


Hd 'V' V 

when Hj > dvr^a^M^, 


H, 


tan tan 


-1 


-27raMo 




(27) 

— Ana'^Mg 

(28) 


where, G = ^dvr^a^M^ — i7|. The presence of tangent function in Eq.(28) implies the peri¬ 
odic variation in 0 when the magnitude of driving held is above the value of 27raMs(=181.6 
Oe) and it is referred to the Walker limit of the held^®. The velocity(Eq.(22)) and the 
width(Eq.(23)) of the domain wall can be obtained for the driving helds corresponding to 
below and above the Walker limit of held from Eqs.(27) and (28) respectively. From Eq.(27), 
we can show that when t oo, the 0(t) reaches the constant value, which can be called as 
saturated excitation angle(0s) and it can derived as 


0s = 0(cxo) = tan 


-1 


2'KaMs 

Hh 



= constant. 


(29) 


Correspondingly, the saturated velocity(ns = n(oo)) and the saturated width(hFs = hF(oo)) 
can be calculated by substituting Eq.(29) in Eqs.(22) and (23) respectively. From Eq.(29) 
we can observe that at t = oo, dp/dt = 0. 

The velocity of the domain wall can also be obtained by reducing the Eq.(4b) at x = X 
and using Eqs.(20) as 


v{t) = 


dX 

dt 


W{t) I Hu + Hy \ dp 

a p Ht-Hy dt 



(30) 


By substituting dp/dt = 0 for f = oo in Eq.(30), we get the saturated velocity of the domain 
wall in the absence of current and transverse magnetic held 


Vs = v{oo) 


iHdWs 

a 


(31) 


For the understanding of domain wall motion, the quantities v{t), W{t)/W{t)) and pit) are 
plotted in the hgures 4(a-c) for the driving helds below the Walker limit(181.6 Oe) and 


16 





























R. Arun, P. Sabareesan and M. Daniel 


the corresponding inset figures exhibit the oscillatory behaviour of the domain wall for the 
driving field above Walker limit. 

From FIG.4(a), it is observed that while applying the driving field along the positive 
x-direction, the domain wall starts to move slowly in the same direction. Also, the initial 
velocity n(0) is consistent with the expression n(0) = a'-^WoHd!{)- + ~ 0 derived from 

Eq.(22). After a time around t=0.5 ns, the velocity of the domain wall gets saturated when 
the reduction in Zeeman energy of the magnetic strip equal to the dissipation of energy 
through damping^®, and the corresponding saturated value varies depending on the strength 
of the driving field which is also consistent with the equation (31). The plots corresponding 
to Hd=50 Oe and -50 Oe imply that the direction of the velocity of the wall is same as the 
direction of the driving field and the magnitude of the velocity is independent of the direction 
of the driving field. The saturated velocity increases with the increase of driving field when 
its strength is lower and decreases with the increase of driving field when the strength of 
driving field reaches the Walker limit is shown in FIG.4(a). This can be understood from the 
expression of the saturated velocity which is proportional to the product of driving field, and 
saturated width and the saturated width decreases with the increase of driving field as shown 
in the FIG.4(b). Initially, the width ratio assumes to be unity and decreases with respect 
to time and reaches its saturation Ws/W{Q) at around t=0.5 ns for different driving fields. 
The time variation of the width ratio is irrespective of the direction of the driving field. 
The saturated width ratio 144/14^(0) decreases with increase of driving field implies that the 
saturated width Wg decreases when the driving field increases as shown in FIG.4(b). The 
excitation angle (j) starts from zero and reaches its saturation (pg around time t=0.5 ns for 
different driving field is shown in FIG.4(c). The variation in the sign of the excitation angle 
represents that the the magnetic moments in the domain wall excite upwards or downwards 
with respect to the plane of the nanostrip corresponding to positive or negative direction 
of the driving field respectively. The decrease in the width ratio(see FIG.4(b)) and increase 
in the excitation angle((see FIG.4(c)) with respect to time implicitly indicate the increase 
in distortion observed in the domain wall for all driving fields below the Walker limit. The 
distortion of the domain wall is associated with damping and it increases with time upto 
around 0.5 ns after that the rate of decrease in zeeman energy balances the rate of damping 
in energy dissipation and the wall starts to move with constant velocity. 
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2. In the presence of transverse magnetic field alone(Hy 7 ^ 0 and b = = 0): 

The variation of the excitation angle of the domain wall after applying the transverse 
magnetic held is obtained by substituting b = Hd = 0 in Eq.(21) is given as 


(1 + = —a7(27rMs sin 20 + Hysmcl)). 


dt 


Eq.(32) can be integrated as 


dc() 


ay 


4 ( 0 ) 


2 'kMs sin 2 (j) + Hy sin 0 


1 + a^ 


-dt. 


We obtain the following equation after integrating the Eq.(33). 


tan ■ 




Hy + AttMs cos 0 
sin 0 


AnMs 


= exp In 

exp ( Hy\n tan 


Hy + 47 rMs cos 0 ( 0 ) 


sin 0 ( 0 ) 


(32) 


(33) 


0 ( 0 ) 


exp 


a7(Tr2 levr^M^) 


1 + a^ 


(34) 


By applying the initial condition 0(0) = 0 in the above equation (Eq.(34)), 0(t) is obtained 
as 


0(t) = 0. 


(35) 


On substituting Eq.(35) in Eqs.(22) and (23), we get 


v{t) = 0 and W{t) = Wo/\ll - ( ^ 


= constant. 


(36) 


Eqs.(35) and (36) represent that there is no variation in the excitation angle, velocity and 
width of the domain wall with respect to time only when a transverse magnetic held is 
applied. 


E. Analytical solution for the domain wall motion using small angle 
approximation 

The analytical solution for the excitation angle 0 under the small angle approximation are 
valid only for below the Walker limit. However for above the Walker limit, the domain wall 
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would exhibits oscillatory behaviour which has been observed from Eq.(28) in the previous 
section. In the forthcominng section, we discuss the analytical solution for the excitation 
angle, the velocity and the width of the domain wall in presence of transverse magnetic held 
along with driving held and current. 


1. The general case(b, Hy 0) 


To understand the motion of the domain wall driven by the driving held and current in 
the presence of transverse magnetic held, the excitation angle of the wall is obtained by 
rewriting the Eq.(21) using the small angle approximations sin0 0 and cos0 ps 1 as 


(1 + = 'yHd - a'j{4:TrMs + Hy)(j) + 


hPn 


Hk 


1 + 


47rM,02 _ 

K 


1/2 


(37) 

Since, 0 is very small, Eq.(37) is rewritten after making the approximation [1 + — Hy)/ 

[1 + (47rM,02 _ Hy)/2Hk\ as 


^ = E 4,{tY - F m + G, 


(38) 


where E, F and G are given by 

(39a) 
(39b) 
(39c) 

Eq.(38) is in the form of the well known Riccati equation, the solution can be written as 

where (l)p is a particular solution of Eq. (38) and Z(t) is yet to be determined. Since the time 
independent solution of Eq.(38) is a particular solution of the same equation, 0p is derived 


E = 


F = 


27rM,(«-0^^1-5^ 

HkWo{l + « 2 ) ^ 

Q!7(47rMs + Hy) 


1 + 


G = 


1 + ci^ 


iHd 


[a 


Ob 




Hk 


Hy 

2H, 
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after substituting 0=0p=constant in Eq.(38). The result reads 
^ ^ + Hy)Wo 


AirMsia - i)h^l - ^ 


1 ± 


SirMsia — Ci^\ 1 — 


Hy 

Hk 


aV^fc(47rM, + i/p)Wo 


iH, 


IRn 


Hk 


Hy 


2Hu 


(41) 


Where 0+ and are two time independent solutions of Eq.(38). On substituting Eq.(40) 
in Eq.(38), Z{t) is derived as 

dZ{t) 


dt 


= {F- 2E(f)p)Z{t) - E. 


(42) 


Integrating Eq.(42), we get 

Z(() = 


+ Csexp [(F - 2E4tp)t ], 


F — 2Ecf)p 

where Os is the constant of integration. After substituting Eq.(43) in Eq.(40), we obtain 

E 


(43) 


0(t) = 


+ O 3 exp [(F - 2 F 0 p)t] 


-1 


+ 


xp. 


(44) 


F — 2E(j)p 

The constant of integration O 3 is evaluated using the initial condition t = 0, </> = (/)o in 
Eq.(44). The result reads 

^ ^ ^ (45) 


\F-2E(t)p (j)o-y.p^ 

By substituting the equations (45), F = E{(j)^ + 0_) and 0p = 0+ in Eq.(44), we get 

(00 - 0 +)( 0 - - 0 +) 


0 (^) = S 0 + + 


, 0 _ < 0 +, & 7 ^ 0 . 


(46) 


(00 - 0+) - (00 - 0-) exp(F(0_ - 0+)t) 

Eq.(46) gives the analytical solution for the excitation angle of the domain wall under small 
angle approximation in the presence of a transverse magnetic field{Hy), driving field(i7c;) 
and spin-transfer torque(&). The velocity and width of the domain wall can be determined 
from Eq.(46) using Eqs.(22) and (23) respectively. Since, F(0_ — 0_|_) < 0, one can show 
from Eq.(46) that the excitation angle 0 reaches the saturated value 0^ when t —)■ 00 . 


0 (cxo) = 0s = = constant. 


(47) 
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Eq.(47) implies that the excitation angle saturates as time passes and simillarly the veloc- 
ity(Eq.(22)) and the width(Eq.(23)) of the domain wall reach the saturated values Vs and 
Ws respectively because they are explicitly depend on 0. By substituting d(j)/dt=0 and 
W{t) = Ws in Eq.(30), we can derive the saturated velocity as 


_ 7HdWs Hk + Hy f 

a \j Hk- Hy a 

From the saturated velocity expression, we observed that the saturated velocity of a domain 
wall is high only when the current and the driving held are applied in the opposite directions 
whereas it is low for the current and the driving held are applied in the same directions. In 
the absence of driving held, the transverse magnetic held does not inhuence the saturated 
velocity of the domain wall and it is controlled by the adiabatic spin-transfer torque b. In the 
presence of driving held, the saturated velocity increases with the transverse magnetic held. 
It implies the saturated velocity of a domain wall cannot be increased by the transverse 
magnetic held in the absence of driving held. 



2. For the case when Hd, Hy ^ 0 and 6 = 0 


When the current is switched oh(6 = 0), the excitation angle (p given in Eq.(46) is 
undetermined. To hnd the solution of the excitation angle in the presence of transverse 
magnetic held, driving held and absence of current, Eq.(38) is written after substituting 
6 = 0 as 


dp _ -fHd a-f{47iMs + Hy) 

dt~{l + a^) (l + a2) 


On integrating Eq.(49), we get 

In {Hd - «(47rM, + Hy)P) = + Hy)t ^ 

(1 + W) 

where, the constant of integration 04 is obtained using the initial condition t 


(49) 


(50) 


0 , •)> = '(> 0 . 


C 4 ^ hi - a( 47 rA/s + i/„) 0 o) ■ 


(51) 


After substituting C, in Eq.(50), the value of (p for the domain wall driven by driving field 

in the presence of transverse magnetic held is derived as follows 


p{t) = 


Hd 


a{ATTMs + Hy 


+ 00 - 


Hd 


a{ATrMs + Hy 


exp 


Q;7(47rMs + Hy)t 

1 -f 


(52) 
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The corresponding velocity and width of the domain wall can be found from Eqs.(22) and 
(23) respectively using Eq.(52). The saturated excitation angle(0s) is given from Eq.(52) at 
t = oo a.s 


a{47rMs + Hy)' 

The Eq.(53) shows that the increase in transverse magnetic field decreases the magnitude 
of saturated excitation angle. This implies that the excitation of the domain wall from the 
plane of strip can be controlled by the transverse magnetic field. 


F. Displacement of the domain wall 

The displacement of the domain wall X{t) is obtained by integrating Eq.(30) with respect 
to time and using the initial condition 0(0) = 0,X(0) = 0. 


^ f(54) 


Eq.(54) can be rewritten using Eq.(23) as 


X{t) = 




{iH, - f) dt 


a 


T _ Ejl Jo 

Hk 


l + ^sin^0-^ 


Hk 


Hk 


COS( 


-6- / dt. 
a Jo 


(55) 


As it is difficult to integrate the Eq.(55) in its present form, we implement the small angle 
approximations sin 0 ps 0 and cos 0 1 and the result reads 


X{t) = 






{-fHdh - h) - 6- / dt. 


a 


(56) 


Where, the integrals R and I 2 are given by. 


h = 


dt 


h = 


d(j) 




(57) 
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The integrals Ji and I 2 are evaluated using Mathematica 8.0 and the results read 


h = 




EFiF2{(j)+ - (j)-){ 


IfFi 


(0+ - 


+ Fi log [F^\{(p- - <p+Fs){Hk - Hy) - 47rM,(F3 - 1)0'_0+ + F 2 F^)\ 

+ Fi log [F^\{(t>- - ct>+F^){Hy - Hk) + 47rM,(F3 - 1)0+0- + F^F^)] 

(58a) 


F = 


Hk 


AtiM, 


log 


H ^ 
Hk, 


AnMc: 




+ w 1 + 


47rMo 


H,. 


-d>2_ 

Hk ^ Hk 


(58b) 


where, 


Fi = ^Hk-Hy + 47rM,0^, 


(59a) 

F 2 = ^Hk-Hy + 47rM,02_, ( 59 b) 

F 3 = exp {F(0+ - 0_)t} , (59c) 

F 4 = ^(Ffc - F,)(0_ - 0+F3)2 + 47rM,(F3 - 1)20^02_ (59d) 

and F is given by (39a). By substituting the Eqs.(58) in Eq.(56) gives the expression for 

the displacement of the domain wall. The above result is valid only for bT^O, if b=0, the 
integrals R and I 2 are undetermined. Hence, to calculate the displacement of the domain 
wall for the case of b=0, the integrals R and I 2 in Eq.(56) are evaluated using the same 
initial condition 0(0) = 0,X(0) = 0 with the excitation angle given by Eq.(52) instead of 
Eq.(46). Therefore, the integral Ji is given by 

Ji =^—^^a'y{A7rMs + Hy)t + (1 + a^) 

log [G'r'( 47 rM,(Gi - 1)H^^ + aGi{Hy + A7rM,)[a{Hy - Hk){Hy + dvrM,) - /^G' 2 G 3 ]) 

(60) 

and the integral I 2 is same as given in Eq.(58b) except that 0 is obtained from Eq.(52). 


where 


Gi = exp 


a'y^Hy + A7iMs)t 


1 + a^ 

G2 = + {Hk - Hy){ATiMs + HyYa\ 

V4(Fi - IYtiM.HI - Ffi^Hy - Hk){Hy + 47rM,)2«2 


G.= 


aF^^k{Hy + 47rM,) 


(61) 
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IV. CONFIRMATION OF ANALYTICAL SOLUTIONS WITH 
NUMERICAL RESULTS AND DISCUSSIONS 

In this section, the excitation angle, velocity, width ratio and displacement of the domain 
wall are obtained from the approximated analytical solutions Eqs.(46) and (52) are verified 
with the corresponding numerical results obtained by integrating the dynamical equation 
(21) using Runge-Kutta-4 algorithm with the initial condition 0(0) = 0 and by using the 
experimentally measured values of the material parameters of Cobalt nanostripes as given 
by Ms = 14.46 x 10^ Am-0 47 rM, = 1.8 x 10^ Oe, A = 2 x lO'^^ Jm-0 7 = 1.9 x 
10^ Oe-^s-\ Hk = 500 Oe, P = 0.35 and ^ = O.OP'^^’^^ 


Case A: Effect of the current in the Domain Wall motion(6, c 7 ^ 0; Hd = Hy = 0) 

In order to substantiate the statement that the moving domain wall maintains the con¬ 
stant velocity and width, the analytical and numerical solutions for the velocity(n), width 
ratio(lU/lU(0)), excitation angle(0) and displacement (A) of the domain wall are plotted 
against time in FIGs.5(a-d) respectively for the different values of current(&=-200,-400,- 
600,-800 and 200 m/s) in the absence of driving field and transverse magnetic held. The 
open circle represents the numerical results and solid line represents the analytical results ob¬ 
tained from Eq.(21) and Eq.(46) respectively. The quantities 0,n and 1U/1U(0) are initially 
changing with time and after a fraction of nanosecond they reach the saturated values such 
as 05, Us, 1E//IU(0) as shown in FIGs.5(a-c). This is because of the equivalence between the 
rate of energy supplied by the incoming electrons and the rate of energy dissipation through 
damping. 

In FIG.5(a) we observe, initially the velocity of the domain wall seems to be equivalent to 
—b, and this can be explained by calculating the initial velocity n( 0 ) = —&(l-|-a^)/(l-|-a^) ~ 
—b from Eq.(22) using the initial condition 0(0) = 0. When the driving held and the 
transverse magnetic held are switched oh, the saturated velocity of the domain wall is given 
by Vg = —c/a (Eq.(31)). Hence, the adiabatic spin-transfer torque is considered to be most 
important for the initial velocity of the domain wall, whereas, the nonadiabatic spin-transfer 
torque controls the hnal velocity of the domain wall. 

Eventhough, the initial velocity is maximum and completely depends on the adiabatic 
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FIG. 5. (color online). The time variation of (a) velocity v{t), (b) width W{t), (c) excitation angle 
(p{t) and (d) displacement X{t) of the Cobalt domain wall for different values of the current(6 =- 
200,-400,-600,-800 and 200 m/s) without driving field and transverse magnetic field. The same 
colour of the solid line and open circle represent the analytical and numerical results respectively 
for the given spin-transfer torque. 


spin-transfer torque(6), b itself cannot move the domain wall for a long distance, because 
the saturated velocity becomes zero without nonadiabatic spin-transfer torque. The reason 
behind the vanishing of saturated velocity in the absence of nonadiabatic spin-transfer torque 
is the absorbance of the adiabatic spin angular momentum of the incoming electrons by the 
distortion of the domain wall, caused by the adiabatic spin-transfer torque, so that the net 
adiabatic torque on the wall vanishes and as a consequence the wall stops. The nonadiabatic 
spin-transfer toruqe(c) is responsible for the nonzero final velocity of the domain wall, which 
behaves as a nonuniform magnetic held that can sustain a steady state motion of the wall. 
Eventhough the magnitude of c is about 2 orders smaller than the magnitude of b, the 
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FIG. 6. (color online), (a) Spatial profile of the normalized x-component of magnetization M^/Mg 
of the wall for different time, (b) Velocity of the wall with respect to wall position and (c) Normal¬ 
ized z-component of the magnetization of the domain wall at the center of the wall as a function of 
wall position, for the value of adibatic spin-transfer torque b = —800m/s along with nonadiabatic 
spin transfer torque. 


saturated velocity is large because it is inversely proportional to the damping parameter 
which is very small. The reason for the initial velocity to be larger than the saturated 
velocity is due to the fact that the spin current supplies energy to the domain wall which 
is greater than the damping of the wall at the initial time. As time increases, the damping 
increases through the distortion of the wall width and excitation of the domain wall. When 
the supplied energy equals the damping energy, the wall moves at a constant velocity which 
is less than initial velocity"^®. It may be noted that, the direction of the velocity of the wall 
is always opposite to the direction of current, which means that the domain wall is dragged 
towards the direction of flow of electrons. Similar to the velocity, the width ratio of the 
domain wall also decreases with time for all values of current below the Walker limit and 
reaches the saturated value Wg after a fraction of nanosecond. However, when the current 
increases, the saturated value of width ratio decreases and the variation of the width is 
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independent of the direction of the current as shown FIG.5(b). 

The plot for the excitation angle (j){t) indicates that while the domain wall moves, the 
magnetic moments pointed in the positive x-direction precesses from -x to +x direction 
in a plane which is inclined by an angle (j) with positive y-direction is shown in FIG.5(c). 
The change in the excitation angle indicates that the plane in which the precession of the 
magnetic moments of the domain wall takes place due to the adiabatic and nonadiabatic 
spin-transfer torques is forced out of the plane of the strip when the current is applied. The 
excitation angle increases with time and reaches the saturated value (ps, which also increases 
with the current. The saturation is attained due to the demagnetization field which tries 
to pull the magnetic moments from out of plane to inplane. When b > 0 and & < 0 the 
magnetic moments of domain wall excite in the upward and downward direction to the 
plane of the strip respectively. Using the numerical integration of Eq. (22) (open circle) and 
analytical solution of Eqs. (54)-(59) (solid line), the displacement of the domain wall is plotted 
in FIG.5(d). Initially, the displacement of the domain wall is not linear with respect to time 
since the velocity is not a constant and later, the displacement becomes linear because of 
the velocity is constant. And also it is observed that the domain wall is displaced in the 
direction opposite to the direction of the current. From the plots corresponding to b =200 
and -200 m/ in FIGs.5(a),(c) and (d), we observed that when the direction of current is 
reversed, the dynamical quantities velocity, the excitation angle and the displacement X 
exhibit the same variation with time however with a reversed sign. 

The spatial variation of the normalized x-component of magnetization M^/Ms at different 
times(t=0.0,0.5,1.0,2.0,3.0 ns) against the position of the center of the wall, to observe the 
displacement, velocity and distortion of the domain wall have been plotted in FIGs.6(a-c) 
respectively. From FIG.6(a), it is observed that the wall is not stopped and moves with a 
constant velocity for a long distance due the nonadiabatic spin-transfer torque. At x = 0, the 
velocity of the wall is maximum and then decreases untill it reaches the saturated velocity 
and thereafter it remains constant as shown in FIG.6(b). While the decrease in the velocity 
of the wall is due to the adiabatic spin angular momentum, its constancy is due to the 
nonadiabatic spin angular momentum^^. In FIG.6(c), the z-component of the normalized 
magnetization(M 2 /Ms) at the center of the domain wall has been plotted as a function of 
the position. The M^/Ms implies the excitation of the domain wall from the plane of the 
strip and it shows that the excitation of the domain wall increases initially and maintains 


27 


R. Arun, P. Sabareesan and M. Daniel 


as constant for a long distance. Moreover, the saturation of excitation is very small, about 
\Mz\/Ms K, 0.06, and one can understand that the domain wall is still of Neel wall type with 
small excitation. 


Case B: Effect of the transverse magnetic field on current and field driven 
domain wall motion(6, c, Hy ^ 0) 



FIG. 7. (color online). The variation of the saturated velocity against the transverse magnetic 
field for different strengths of current and driving field in different directions have been plotted. 
The solid lines have been used for analytical results and the corresponding circles of the same color 
have been used for numerical results. 

From the previous sections, it is observed that the initial velocity of the domain wall 
is controlled by the adiabatic spin-transfer torque and the hnal velocity is controlled by 
nonadiabatic spin-transfer torque and driving field. Further, we observed that the domain 
wall does not move and excite when the transverse magnetic held alone is applied. In this 
section, we study the impact of the transverse magnetic held on the saturated velocity of 
the domain wall driven by driving held and current and the results are plotted in FIG.7 
between saturated velocity and transverse magnetic held for the diherent values of current 
and driving held. The solid line represents the analytical results which is plotted by using 
Eqs.(46) and (48) whereas the open circle indicates the numerical results which is obtained 
by solving Eqs.(21) and (30) numerically. 

The plots corresponding to Hd = 50 Oe, 6 = 0 m/s and Hd = —50 Oe, 6 = 0 m/s show 
that the saturated velocity increases from 857 m/s to 2666 m/s and -857 m/s to -2666 m/s 
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respectively when the transverse magnetic held is increased from 0 Oe to 400 Oe. Here, 
the saturated velocity is increased by the transverse magnetic held, because the magnetic 
energy of the domain wall decreases by the transverse magnetic held and the wall moves 
faster in order to have a larger rate of energy damping to dissipate the energy fastly."^^’^® 
The saturated velocity corresponding to the plots oi b = —800 m/s, Hd = 50 Oe and 
b = 800 m/s, Hd = —50 Oe increases from 1400 m/s to 3374 m/s and -1400 m/s to - 
3374 m/s respectively when the transverse magnetic held is increased from 0 Oe to 400 Oe. 
Similarly, the saturated velocity corresponding to the plots b = 800 m/s, Hd = 50 Oe and 
b = —800 m/s, Hd = —50 Oe increases from 259 m/s to 1944 m/s and -259 m/s to -1944 m/s 
respectively when the transverse magnetic held is increased from 0 Oe to 400 Oe. These plots 
imply that the saturated velocity can be increased by the transverse magnetic held and it is 
irrespective of the directions of the driving held and current. However, the saturated velocity 
is higher (lower) when the directions of driving held and current are antiparallel(parallel). 
The reason for the suppression of saturated velocity when current and driving held are 
applied in the same direction, can be explained as follows. In general, a domain wall moves 
in the opposite direction of current and same direction of driving held. Therefore, when both 
the current and the held are applied in the same direction, the motion by the current (driving 
held) will be opposed by the driving held(current) so that the velocity of the domain wall 
decreases. 

In the above two cases(A and B), from the hgures 5(a-d), 6(a-c) and 7, we can observe 
that the plots corresponding to the solid line resembles with the plots corresponding to the 
open circle. This indicates that the obtained approximated analytical solutions given in 
Eqs.(46) and (52) match with the numerical results. Therefore, these analytical solutions 
can be used to understand the dynamics of domain wall in ferromagnetic nanostrip for the 
diherent materials. 


V. CONCLUSIONS 

In the present paper, the dynamics of transverse Neel wall in a ferromagnetic nanostrip 
in the presence of current, driving held and transverse magnetic held is studied by solving 
the Landau-Lifshitz-Gilbert equation with spin-transfer torques. By converting the LLG 
equation into spherical coordinates and using trial functions, the domain wall parameters 
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such as the excitation angle, velocity, width and displacement are obtained. Under small 
angle approximation, the equation for cj) is further reduced to the form of the Riccati equation 
and using its solution, the above mentioned four quantities in the presence of current, driving 
field and transverse magnetic field are obtained analytically. The results show that the initial 
velocity of the domain wall can be controlled by the adiabatic spin-transfer torque and the 
saturated velocity can be controlled by the nonadiabatic spin-transfer torque and driving 
field. The direction of the saturated velocity is antiparallel and parallel to the directions of 
the current and the driving held respectively. For the current driven domain wall motion 
with nonadiabatic spin-transfer torque, the saturated velocity of the wall linearly increases 
with the current. But in the case of held driven domain wall motion, the saturated velocity 
no longer increases with the driving held, because the saturated velocity is proportional 
to the product of the driving held and the saturated width decreases with increase in the 
driving held. In the presence of a transverse magnetic held alone the domain wall is at rest, 
but the width of the wall is increased at the time of held which is applied and thereafter the 
width is also constant. 

Numerical results showed that the domain wall is driven by the current in the presence 
of the transverse magnetic held, the saturated velocity remains constant. However the 
initial velocity can be increased with transverse magnetic held when the initial value of the 
excitation angle 0 is slightly perturbed. Whereas for the held driven domain wall motion, 
the saturated velocity of the wall increases considerably, when the transverse magnetic 
held is increased. The transverse magnetic held is increased from 0 Oe to 400 Oe, the 
corresponding saturated velocity increases from 857 m/s to 2666 m/s, for the case of = 50 
Oe, 6 = 0 m/s and from 1400 m/s to 3374 m/s, for the case of = 50 Oe, b = —800 
m/s. Further, the numerical results showed that the saturated velocity is increased by the 
transverse magnetic held with the irrespective of the directions of the driving held and 
current and the saturated velocity is higher(lower) when the directions of driving held and 
current are antiparallel(parallel). The obtained approximated analytical solutions given in 
Eqs.(46) and (52) match with the computed numerical results. 

In conclusion, the obtained analytical results of the dynamical parameters namely the 
excitation angle, velocity, width and displacement closely coincide with the numerical results. 
While the transverse magnetic held has no ehect on the saturated velocity of the domain 
wall in the current driven case, in the held driven case the saturated velocity of the domain 
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wall increases by a large amount due to transverse magnetic field. The compatibility of 
analytical results with numerical results implies that the analytical model with small angle 
approximation can be considered for transverse Neel domain wall dynamics in ferromagnetic 
nanostrips. From the above, it is inferred that transverse magnetic field plays a crucial role 
in enhancing the saturated velocity of the domain wall which will be useful to design high 
and low speed domain walls depending upon the real time applications and also promises to 
make the efficient storage devices. 
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